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Abstract. How do charge and density fluctuations compete in ionic fluids near 
gas-liquid criticality when quantum mechanical effects play a role ? To gain some 
insight, long-range <f>^ ± /r d+CT interactions (with a > 0), that encompass van 
der Waals forces (when a = d = 3), have been incorporated in exactly soluble, d- 
dimensional 1:1 ionic spherical models with charges ±<?o and hard-core repulsions. 
In accord with previous work, when d > min{<7, 2} (and qo is not too large), the 
Coulomb interactions do not alter the (qo = 0) critical universality class that 
is characterized by density correlations at criticality decaying as l/r d ~ 2+v with 
ri = max{0, 2 — a}. But screening is now algebraic, the charge-charge correlations 
decaying, in general, only as l/r d+CT+4 ; thus a = 3 faithfully mimics known 
noncritical d = 3 quantal effects. But in the absence of full (+, — ) ion symmetry, 
density and charge fluctuations mix via a transparent mechanism: then the 
screening at criticality is weaker by a factor r 4 ~ 2v . Furthermore, the otherwise 
valid Stillinger-Lovett sum rule fails at criticality whenever rj = (as, e.g., when 
a > 2) although it remains valid if r\ > (as for a < 2 or in real d < 3 Ising-type 
systems). 



An ionic fluid, such as an electrolyte or a plasma, is characterized in thermal 
equilibrium by the screening of the long-range Coulomb interaction potential, 
z T z v qo/r d ~ 2 (say, in d > 2 dimensions) between ions of charges z T qo and z v qo. 
Following Debye-Huckel theory pQ one expects the charge-charge correlation function,! 
G Z z(r), to decay exponentially as e~ r /'> Z oc where, as the overall ion density p 
becomes small, the screening length £z,oc(T,p) should approach the Debye length 
£ D = 1/k d oc ^T/q 2 p. More generally, the screening of an external charge in a 



conductor should be characterized by the Stillinger-Lovett sum rule p] Ej which does 
not require exponential screening [3]; this states that when k — > 0, the charge structure 
factor, essentially the Fourier transform! °f G zz (r), behaves as 



s zz (k) = o+e Ztl k 2 -Q 2 k 4 -ai^ + ... with ui = &>. (i) 



Note also that the vanishing leading term simply reflects the requirement of 
electroneutrality; if the screening is exponential only the further powers k 21 with 
I = 1,2, ... can arise; but, in general, which specific, nonanalytic powers appear and 
with what amplitudes is a matter of prime interest. 

J Our notations, which are fairly standard [l] are set out explicitly in which will be denoted AF. 
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The exponential screening of charge has been proven rigorously at low densities 
for classical systems 0; but sufficiently strong short-distance repulsions between 
oppositely charged ions are essential while any further non-Coulomb interaction 
potentials, say ip TV (r) between ions of species r and v, must be of short-range, 
decaying, e.g., exponentially fast when r — > oo. But then, in allowing for quantum 
mechanics, one should first recognize that real polarizable ions also interact via 
fluctuating induced dipole-dipole or van der Waals forces that fall off only algebraically. 
Specifically, if for generality we consider integrable long-range potentials that decay as 
^^ v /r d+rT (so requiring a > 0) 0, van der Waals forces may be described by d = <r = 3. 
Then one may ask what will be the consequences of such long-range power-law forces 
for exponential screening and for the expansion (JJ. 

However, even in a point-ion model of a plasma, the quantum-mechanical 
fluctuations of position lead to effective particle-particle potentials, as manifest in 
the corresponding correlation functions, that decrease as 1/r 6 (for d = 3) 0|B], just 
as for polarizable ions (or neutral species). Furthermore, it is known that in quantum 
plasmas at low density, the screening of charge is no longer exponential: indeed, G zz (r) 
has recently been shown to decay as 1/r 10 0|Sj. 

On approach to gas-liquid or liquid-liquid criticality, in an ionic system, the 
situation is further complicated because the density fluctuations become divergent and 
it is natural to ask if this should not seriously affect charge screening, the Stillinger- 
Lovett rule, etc. To be more explicit, in a short-range, nonionic system, the density- 
density correlation function G NN (r) decays as e~ r ^ N -°° where, on the critical isochore 
p = Pc, the density correlation length (T, p) diverges as l/t u with t = (T — T c )/T c 
and v > i; but at criticality, one has 

G c NN (r)^D% N /r d - 2+r > when r - oo , (2) 

with D C NN finite and rj > 0. Might not the divergence of £ N ,oo couple in some way to 
the screening length £ Zj00 in a corresponding ($^ = 0) classical plasma and cause it 
to diverge at criticality ? Likewise, might not the slow decay of G NN (r) at criticality 
change the 1/r 10 quantal decay of G zz (r) that has been established at low densities ? 

Of course, the values of the critical exponents v and rj depend on the critical 
universality class, and the influence of Coulomb couplings on criticality has been 
an important experimental and theoretical question ever since seemingly convincing 
observations on certain electrolytes suggested that ionic criticality might realize a 
new or different type of critical behavior [9-11]. While the experimental issues for 
electrolytes may now be regarded as largely settled — in favor of no change in critical 
character — the theoretical situation remains open. And neither experimentally nor 
theoretically has a clear picture of the charge correlations near criticality yet emerged. 

One question of particular relevance fTTJ concerns the role of ion- symmetry. Thus, 
the simplest theoretical description of a 1:1 electrolyte is afforded by the so-called 
restricted primitive model or RPM, in which equisized hard spheres carry charges +qo 
and —go- The RPM is precisely ion symmetric so that it is plausible that charge and 
density fluctuations will remain effectively independent even at criticality. Indeed, 
careful simulations have recently established that criticality in the d — 3 RPM 

still exhibits Ising behavior as do simple, nonionic fluids. But in reality, + and — ions 
are never identical so that ion nonsymmetric models are of especial interest. While 
simulations of 1:1, 2:1 and 3:1 hard sphere models with unequal diameters have been 
undertaken [2], their critical behavior has not, as yet, been resolved. 
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To address some of these problems, we have recently analyzed j2] a class of exactly 
soluble d-dimensional multicomponent spherical models in which particles of different 
species r, v, . . . reside on distinct but equivalent interlaced sublattices with nearest 
neighbor spacing a. Then, in a 1:1 'lattice electrolyte,' nearby ions of species + and 
— with charges ±go avoid collapse because opposite charges are never closer than 
the minimal, intersublattice distance, say ao < a 2\. In addition to the ±q 2 /r d ~ 2 
Coulomb potentials, short-range attractive potentials ip±±(r) of magnitude, say, k B To 
are then introduced; on setting <?o = 0, these suffice to yield standard spherical model 
critical behavior, with rj = etc., at T c (qo = 0) ~ To [2 . 

Two main conclusions emerge from AF 0. First, the character of criticality 
remains unchanged when the charges are switched on (provided qo is not too large). 
Second, ion symmetry plays a crucial role: thus in symmetric models, the charge 
screening length £ z ,oo(T, p) remains of order £ D even at (T c , p c ) although it gains 
singular corrections when t — > 0; but, on the contrary, for asymmetric fluids, 
£z,oo diverges on approach to criticality, precisely matching the density correlation 
length £ N ,oo ~ l/t u . Furthermore, the Stillinger-Lovett rule is then violated at 
criticality. Here we extend AF by incorporating additional, long-range &± ± /r d+<T ion- 
ion potentials in order to cast some light on the further role that quantum mechanical 
fluctuations might play. (To avoid technical complications we suppose a ^ 2, 4, • • • .) 
As we report, the findings prove instructive. 

The necessary analysis follows closely the lines set out in AF: accordingly we focus 
on the principal results sketching only the basic technical points. (Further and fuller 
details will be presented elsewhere QS].) Thus Lagrange multipliers (A + , A_) = A(T, p) 
are introduced to satisfy the spherical model constraints (s 2 ) = l for r = +, — , where 
the s T (R) are the usual unbounded scalar spin variables [16-19]. The free energy then 
follows from the eigenvalues, A JV (k; A) and A z (k; A), of the interaction matrix A with 
elements of the form A TV — [X T S TV — tp TV (k)], where the (p TV are proportional to the 
Fourier transforms of the total interaction potentials ip TV (r) . 

Now, as in AF, the crucial result is that the qo/k 2 Coulomb divergence cancels 
out identically from the first eigenvalue which, at small wave numbers k = |k|, then 
behaves as 

A„(k; A) = \k B T [A(T, p) + (R N k) 2 + (R%kf + ...], (3) 

where A (= A/jo in the notation of AF) is found to vanish on the critical isochore 
near criticality like i 7 , while the finite, nonzero length R N (X;qo) measures the range 
of the short range forces [see AF (25)]. The relative contribution of the long range 
interactions to the density variation is embodied in the effective range Rf, [H]. From 
this form for A N one finds (for q not too large) that the critical behavior is always of 
spherical model form with exponents (3 — \ and 

V (a) = max{0, 2 - a}, 7 = 1 - a = (2 - r,)v = /^f^ > (4) 

HHEIIlUn] while d > mm{a, 2} is needed for T c > [ and we suppose (R%) a > when 
a < 2 ]. For a > 2 these leading exponents are, as well known, independent of a and 
the same as for short-range forces; but see also (20]. However, new correction terms 
varying as t e ° with 6 a = |2— a\v(a) will dominate in all properties when a is close to 2. 

The long range effects of the Coulomb forces appear only in the second eigenvalue 
which, for small k, varies as 

A z (k;A) = ^|| [l + (R z k) 2 + (Rtk) 2+ ° + ...] , (5) 
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with Sd — 2-K d / 2 /T(d/2). The net contribution of the long range forces is now 
represented by R z while R z (\; k/k) is of order a/y/To ~ a/q where the ionicity 

1 = ql/a d - 2 k B T , (6) 

measures the overall strength of the Coulomb interactions near criticality. [It should be 
noted that the ellipses in (j^J and © include both further singular terms and analytic 
terms of order k 4 , k 6 , ■ ■ ■ .] 

Finally, we may express the charge and density structure factors in the transparent 

form 

S NN (k) _ 1 - S v B(k) 8 v B(k) 
k B T/4pa d ~ A N (k; A) + A z (k;A)' 

S zz (k) _ yg(k) l-S v B(k) U 
k B T/Apa d A„(k;A) A z (k;A) ' 

where the basic symmetry parameter S v (= A^ /jo in the notation of AF) vanishes 
linearly with the deviation of the potentials (p TV from precise (or effective) ion 
symmetry. Clearly the density and charge fluctuations in ion symmetric models are 
completely uncoupled (at least in quadratic order): the critical density fluctuations 
are thus driven solely by A N (k) while charge screening is entirely controlled by A z (k) . 
But, charge and density fluctuations mix as soon as ion symmetry is lost: to what 
degree is determined by 

B(k) = Ak 4 a 4 [5 ¥ + {R'ky + ...] jSlTl , (8) 

where i?£ measures the strength of the asymmetric parts of the long range forces (and 
hence vanishes with 5 V ). A decomposition similar to (0 holds for S NZ : see AF. 

The factor k 4 in © implies that the "intrinsic" charge fluctuations contribute 
only weakly to S NN . Consequently, except at (T c ,p c ), one has, neglecting analytic 
background terms of order k 4 , 

where \t ~ lA 7 is the isothermal compressibility, while the new length scales are 

£ N ,i(T,p) = iVA 1/2 and £„, CT (T, p) = R c J\ 1/a . (10) 

Since A vanishes like t 7 on approaching criticality, the density correlation length may 
be identified as £ N (T, p) = max{^ Ki i,^} in full accord with the exponent values 

Likewise, the critical point decay (EJ with r] = max{0, 2 — cr} is readily verified. 
Away from criticality, however, matters are somewhat more subtle since, in general, no 
correlation function, G TV (r), can decay faster than the associated \/r d+<J power-law 
potentials Indeed, the first nonanalyticity in ^^(k) yields the large-r behavior 
|22| whence we find, for fixed (T, p) ^ (T c ,p c ), 

G NN (r;T,p)*-^(^) 2 ° ' ~ _L when cr < 2 , (11a) 

when cr > 2 . (116) 



r 



r 

2 



r d+a 



The appearance of the factor (£ n ,<tA) reflects the correction-to-scaling exponent 
9 a = |2 — identified above. 

Now we may study the charge correlations and ask, to start with, about the 
noncritical behavior. When ion symmetry pertains, i.e. 5 V = 0, the decomposition 
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shows that S zz depends only on A z ; then (0 leads directly to the expansion Q 
with, furthermore, full confirmation of the Stillinger-Lovett (SL) relation 0E1 (since 
n 2 D = Sdpqo/k B T). But, by virtue of the factor fc 4 in JjjJ, this remains true even when 
ion symmetry is absent. In addition, the coefficient of fc 4 in is given by 

^,2^e D [Rl-boSlaAT,p)], (12) 

with bo = 0{a 2 ~' y ). Evidently, this moment of G zz is unconstrained and if S v ^ it 
will change sign when, driven by the coupling to the density fluctuations embodied in 
@, diverges like 1/i 7 as T — ► T c on the critical isochore: see lllflll . 

Lastly, the leading nonanalytic term in S zz (k) is given, in Q , by ip — 4 + a while 
its amplitude is 

^ a = & - hS v (R^r + bo%t£\ (13) 

where b\ = 0{a 2 ~ a ). Clearly, this k A+a term is present whether or not ion symmetry 
pertains; but if 5 V ^ it may change sign and its magnitude will diverge, like l/t 27 , 
when (T, p) — > (T c ,p c ). More strikingly, however, the presence of this term ensures 
the destruction of exponential screening; rather one finds j22] 

G zz( r;T,p)«^(^) 4+ ^-^_. (14) 

Thus near to or far from criticality, long range forces always undermine the standard 
picture of Debye screening. Nevertheless, algebraic screening remains in the charge 
correlations. Indeed, we may say that the l/r d ~ 2 Coulomb potential is screened by 
a factor l/r 6+CT or, equivalently, that the long range l/r d+a potential is screened — 
owing to the requirements of "local electroneutrality" — by the factor 1/r 4 . This is a 
central result of our analysis and it is gratifying that on setting a = d — i it reproduces 
the 1/r 10 screening previously found in a fully quantum mechanical analysis of point- 
charge plasmas [3 E] ■ 

In contrast to the loss of exponential screening, ion symmetry is paramount at 
criticality. Thus we see from (0) and JJJ that all the results l|12JI - ill4H remain uniformly 
valid when (T,p) — > (T c ,p c ) provided S v = 0, i.e., that ion symmetry is valid: one 
need only note that the amplitudes R 2 and -D z ,o-(£z.?) 4+<T in l!12ll and l!14t remain 
finite (and nonzero) at the critical point although the former will pick up a singular 
t 1 ~ a correction (as in AF: see [15]V 

On the other hand, for nonsymmetric fluids at criticality, the mixing of the charge 
and density fluctuations depends strongly on a or, more specifically, on 77(a). Indeed, 
displaying only the leading singular terms, we find 

S zz (k) - e D ,c k " [1 + tl(R>k) 2 - a + 0(k 2 )] for V = 2-<j > , (15a) 
= ^ck 2 [Ec-SliR^y- 2 +0(k 2 )] forr; = 0<a-2,(156) 

where J%~ g = a 2 (fl £)- ff , R^ 2 = & 2 (R^Y / R 4 N , £ c = 1 + S 2 v ~a 2 /R 2 N > 1 with 

a (go) = ay/2/SdIo ■ see AF (29). One sees immediately that the Stillinger-Lovett 
sum rule remains valid whenever a < 2 or r\ > 0. This contrasts with the results 
of AF where, in the absence of long-range forces, one has 77 = and the sum-rule is 
violated, precisely the situation that prevails here when a > 2; then the critical system 
may be regarded as an insulator or, at least, as an anomalous conductor! 
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Table 1. Long-distance behavior of charge correlations Gzz(r) at fixed (T, p). 
The density correlation length £n diverges on the critical isochore p = p c as l/t v 
when t = (T — T c )/T c — > 0. The density correlations decay as l/r d + CT away from 
criticality, but at criticality, since r) = max{0, 2 — <r}, they decrease more slowly 
as l/r d ~ a when a < 2 and l/r d ~ 2 when ct > 2. 



G zz (r) 


ion symmetric 


nonsymmetric 




5 V =0 




(T,P) 




a < 2 cr > 2 


= {TcPc), 


~ 1/V d+CT+4 


~ l/r d+4_cr ~ l/r* +a 


± (TcPc), 


~ l/r d+cr+4 


~{i»/rf- 2r > /r d +*-° ~ (£ N /r) 4 /r d + a 



From (15) we can now deduce the long-range behavior of the critical-point charge- 
charge correlation function, namely 



B%^%R^ 1 1 



-2+17 r 2+cr 



for 7i = 0<cr-2. (166) 



As displayed, the results show that although G c zz (r) is driven by the critical density 
fluctuations, with G% N ~ 1/r +?J , the charge correlations are screened relatively 
more strongly when a > 2 than for a < 2, when 77 > 0. In light of the failure of the 
SL rule in the former case rather than the latter, this is, perhaps, paradoxical. On the 
other hand, one might equally conclude that the long-range 1/r + CT potentials are not 
screened, even algebraically, when a > 2 (with r\ = 0) whereas for a < 2 the long-range 
forces are actually screened by factors 1 / r 2r * = A ~ 2a . From that perspective the validity 
of the SL rule when r\ > seems more natural |llj . 

Finally, similar conclusions can be drawn about the charge-density structure 
factor: see AF. Away from criticality we find 

S NZ (k) = «5 v A- 1 ^fc2[ 1 _ (| n ak y + . . .] ? (17 ) 

where = ££ )(T - (i?£) CT /<5 V and one should recall from JUJ that |C£ i<r (T,p)| a 
l/X(T,p) diverges like 1/t 1 while, by (jHJ, R^ vanishes with 5 V , In real space the 
long-range decay can be written as 

x n (c \4-2^(<r) 
^r;T,p)«-^(%) (18) 



r 



so that the cross-correlation function is evidently screened by the factor 1/r 2 relative 
to the long range potentials. Alternatively, one may say that the l/r d ~ 2 Coulomb 
interaction is screened by a factor l/r 4+CT . As regards the "van der Waals case" 
d = a = 3, we may mention that the charge density induced by an infinitesimal local 
external charge decays as 1/r 8 in a fully quantal point-charge plasma [3|Hj. 

At the critical point itself, one finds that S% z (0) vanishes identically if 77 > 
(cr < 2) but that, as in the short-range case j2], S% z (0) = 8 V {^ D / R N ) 2 ^ 0, whenever 
77 = (cr > 2). Moreover, when r — > 00 one has G c NZ lr) ~ 5 ip /r d+ \ a ~ 2 \ so that 
the infection of the charge correlations by the density fluctuations again reduces the 
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screening by a factor r 4-2 '' that is greatest when 77 = 0: Compare with the results for 
G zz (r) as displayed in Table 

In conclusion, we have analyzed the interplay between long-range density 
fluctuations, Coulomb interactions, and power-law \/r d+a (e.g., van der Waals) 
interactions, away from, close to and at criticality in ionic fluids, on the basis of 
two-component, d-dimensional spherical models with hard-core interspecies repulsions 
(2j. Throughout the fluid phase, including low-densities, the power-law forces 
destroy the usual Debye exponential screening. However, algebraic charge-charge 
and density-charge screening is still present: explicitly, although the density-density 
correlation function G NN (r) decays no faster than the interaction potentials |21| . i.e., 
as l/r d+<T , the charge-charge correlation function, G zz (r), decays at large distances 
as D z (T,p)/r d+tT+4 . If one sets a = 3 to mimic the 1/r 6 particle-particle interactions 
that arise in e? = 3 dimensions from quantal fluctuations pJIEj, this result is, indeed, 
in accord with exact results for G zz (r) in a fully quantum-mechanical, point-charge 
plasma at low densities |H] . 

In the critical region, the Coulomb interactions leave the universality class of the 
spherical models unchanged since, as in AF, they still cancel out of the fluctuation 
factor that drives criticality. Nevertheless, radical changes arise in the charge- 
charge and charge-density fluctuations, whether the system is ion symmetric or, more 
realistically, nonsymmetric. The behavior is enforced by a general decomposition 
of the structure factors: see Q and [2]. For ion symmetric fluids, the asymptotic 
amplitude D z (T, p) is always finite and the Stillinger-Lovett (SL) sum rule is satisfied 
even at criticality. For asymmetric systems near criticality, D z is driven by the density 
fluctuations and hence diverges as 1/t 2 "' when p — p c while the SL sum rule remains 
valid. At criticality, however, the density correlations, decaying now as 1/r 2+TI , 
weaken the charge screening still further, more strongly when rj = (a > 2) than if 
r\ > (<7 < 2); see the summary in Table I. Finally, the SL sum rule, characteristic of 
conductors, is satisfied at criticality when 77 > 0, but is violated when 77 = 0, a result 
that may well have validity beyond the ionic spherical models studied here JU- A 
further challenge is to see how far similar results might be obtained for intrinsically 
quantal spherical models such as have been advanced in the past [18,23-25]. 
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